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Abstract 

We provide high-order approximations to periodic travelling wave profiles and to the 
velocity field and the pressure beneath the waves, in flows with constant vorticity over a 
flat bed. 

Keywords'. Travelling water waves, vorticity, velocity field, pressure. 

AMS Subject Classifications (2010): 35Q31, 35Q35, 76D33. 

1 Introduction 

Winds in offshore storms transmit energy to the ocean surface, creating waves. Once the 
waves leave the storm area they become organised into two-dimensional regular wave trains, 
in the form of a regular profile that propagates practically unchanged in a fixed direction. 
These periodic travelling waves are called ‘swell’ in oceanography. Their shape is influenced 
by the underlying currents. The most significant currents on areas of the continental shelf 
and in many coastal inlets are the tidal currents 1201 . These alternating horizontal movements 
of water are created by the gravitational pull of the moon, and to a lesser degree, the sun, on 
the earth’s surface, being associated with the rise and fall of the tide: the current associated 
with a rising water level, called the flood, is directed towards the shore, while the current 
associated with a receding water level, called the ebb, is directed back out to sea. Flows of 
constant vorticity with a flat free surface provide adequate descriptions of pure tidal currents, 
cf. the discussion in [ 6 ], positive constant vorticity 7 > 0 being appropriate for the modelling 
of the ebb current and negative constant vorticity 7 < 0 for the flood current. A spectacular 
example of wave-current interaction is the Columbia River entrance, made by appreciable 
tidal currents one of the most hazardous navigational regions in the world since wave heights 
can easily be doubled in just a few hours m- At this location, tidal velocities in excess 
of 2 m/s are encountered and in winter wave heights in excess of 6 m, up to 14-15 m, are 
common, cf. m- Note that moderate tidal currents reach speeds of up to 0.7 m/s, cf. [36] . 
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while tidal currents with speeds of 5.5 m/s are encountered in between the Scottish mainland 
and the Orkney Islands, in the Pentland Firth — see the data provided in m- 

In this paper we provide accurate approximations for the interaction of waves with under¬ 
lying currents of constant vorticity in the absence of flow-reversal; for theoretical studies and 
numerical simulations of wave-current interactions with flow-reversal we refer to OlISlES] 
and 1271E3], respectively. Let us specify that in a flow where both waves and currents are 
present, given velocity measurements at one point, the ‘current’ is defined as the average 
velocity, and the periodic components that vary around this average are ascribed to the 
wave motion. Since irrotational flows do not present flow-reversals as the underlying current 
must be uniform, one of our main purposes will be to highlight the effects of vorticity on 
the surface wave profile. This problem is of great practical relevance since the detection of 
non-uniform underlying currents from the surface wave pattern is particularly important to 
anticipate and avoid hazardous conditions, where possible. In addition to the wave profile, 
we will also provide approximations for the velocity profile and the pressure throughout the 
fluid. Not only is there a strong interplay of these flow characteristics that is very useful in 
qualitative studies — see um, but in practice information on the state of the sea surface is 
often gathered from subsurface pressure and/or velocity measurements, cf. the discussions in 
[SmiEllMlEHllsll. Let us point out the following counter-intuitive fact: periodic travelling 
waves that propagate at the surface of water with a flat bed in a flow of constant vorticity 
must be symmetric if no flow-reversal occurs and if the wave profile is monotone between 
successive crests and troughs, cf. ilElES]. This means that an underlying non-uniform 
current of constant vorticity does not break the symmetry of irrotational wave trains, so 
that the manifestation of vorticity on the surface wave pattern must take subtler forms. 

The fact that even when both the wave motion and the underlying current of constant 
vorticity are known with accuracy, their interaction produces a significantly different effect 
from that obtained by simply adding the effect of the waves and the currents considered 
separately, cf. the discussions in 1231301131], shows the importance of dealing simultaneously 
with these two flow components. Linear theory (see the discussion in 13) provides the 
dispersion relation 

c = uq — i ^ J 7 ^ tanh^(A;(i) -h 4:gk tanh(A:d) ( 1 . 1 ) 

ZKi Z/C ’ 

for the wave speed c, expressed in terms of the surface current speed i/q in n flow of constant 

27r 

vorticity 7 in water of mean depth d; here = — is the frequency of a wave with wavelength 

L. Note that 0 permits us to understand the effect of vorticity on the propagation speed 
of waves of small amplitude (that are realistically described by linear theory). For example, 
from the point of view of a fixed observer noticing right-propagating waves that interact with 
a current with vanishing surface speed (a setting corresponding to c > 0 and = 0 ), from 
0 we infer that the propagation speed is 

0 ( 7 ) = ^ Jry 2 tanh^(A;d) -h 4:gk tanh(A:d), ( 1 . 2 ) 

irrespective of the sign of 7 . In this setting the underlying current beneath the flat free 
surface y — Q and above the bed y — —d is given by u{y) — 7 ^, so that 7 > 0 corresponds 
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to an adverse current since everywhere beneath the surface the current velocity opposes 
the direction of wave propagation, 7 < 0 corresponds to a favourable current (the current 
velocity points everywhere beneath the surface in the direction of wave propagation), while 
the irrotational case 7 = 0 is characterized by the absence of an underlying current. A simple 
analysis of 0 confirms that 


c(-| 7 |) > c( 0 ) > c(| 7 |) 

for every 7 ^ 0 . This means that favourable currents enhance the wave speed, while adverse 
currents reduce it. In contrast to this analysis of the wave speed, the first-order approx¬ 
imation that is pursued within the framework of linear theory is not conclusive in regard 
to the effect of underlying currents on the wave profiles: at first order, periodic travelling 
waves in flows of constant vorticity are sinusoidal (see i)- Consequently, nonlinear effects 
have to be accounted for and higher-order approximations are needed. We present in this 
paper an approach that provides the second and third order asymptotic expansion of waves 
of small amplitude, expressed in terms of a suitable amplitude parameter b. Following this 
procedure, we can readily obtain other characteristics of the water flow, such as the velocity 
field and the pressure beneath the wave. The knowledge of accurate approximations for 
these wave characteristics allows us to compare the results with those for irrotational flows. 
In the last part of the work we illustrate these flow characteristics for several different types 
of wave-current interactions. Laboratory experiments and numerical simulations for irrota¬ 
tional waves are discussed in ID Elsa, while this type of studies for wave-current interactions 
in flows of constant vorticity were pursued in [HI [161 ESI ESI ESI [30] . The results presented 
below permit a more detailed analysis and numerical simulation. 

2 Mathematical Formulation 

In this paper we discuss two dimensional, periodic travelling waves with constant vorticity. 
We assume that the water is incompressible and inviscid, over a flat bed and acted upon by 
gravity g. In what follows we make no shallowness or small amplitude approximation. In 
order to avoid extended revision we refer to the recent works mini IE] for the mathematical 
formulation of the problem. In this section we make only a brief review of the problem, in 
order to provide sufficient terminology for the fundamental Boundary Value Problem (BVP) 

We are studying two-dimensional waves travelling at constant speed c. This means that 
in a two-dimensional frame moving with the constant speed c, the flow pattern — and, in 
particular, the shape of the surface of the fluid — does not change over time. 

These assumptions, see HD, allow the definitions of the free surface profile by 

S = {(x, y) : —TV < X < TV and y — '^(x)} , 


and the flat bed by 

B = {(x, y) : —TV < X < TV and y — -d) , 
with d > 0 , for the normalized wavelength 2tv. 
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Without loss of generality, the assumption that the waves oscillate around the flat free 
surface ?/ = 0 is made, i.e., 

/ TT 

r]{x)dx = 0. (2.1) 

-TT 

We denote the velocity field of the flow by {u{x^y)^v{x^y))^ with (x^y) G P, where 

V = {(x, y) : —TT < X < TT and — d < y < r]{x)} . 

In this paper we follow the assumption made in that there are no flow-reversals, 
which is formulated as the condition 


u < c throughout the fluid. (2.2) 

We denote the pressure in the fluid by P(x, ^), with (x, y) G V. Neglecting the effects of 
surface tension — a hypothesis that is appropriate for waves of moderate and large amplitude 
(see the discussion in [6]) — we impose that the water pressure is constant on P(x, r]{x)) = 
Patm-) where Patm is the atmospheric pressure. 

Let us recall the following definitions: 

• The vorticity of the flow is defined by 


ry :=Uy- Vx, 

which for the rest of the paper we will assume constant. 
The relative mass flu:xQ is defined by 


(2.3) 


r^{x) 




/ {u{x,y) -c)dy < 0, 

J-d 


(2.4) 


which in fact is independent of x, see m- Moreover, ( |2.4[ ) and ( |2.2| ) show that po 
is negative. This relation shows that the amount of water passing any vertical line is 
constant throughout V. 


• The stream function ^(x^y) is defined as the unique solution of the differential equa¬ 
tions 

= u — c in D, (2-5) 

subject to 

-Ipix, -d) = -PQ . (2.6) 


• In [12] it was proven that the expression 


[u — c)^ + 
2 


+ gy + P -1^ 


(2.7) 


equals a constant E throughout V. The constant Q = E — Patm is called the hydraulic 
head. 

^The terminology ‘relative’ is due to the fact that (u — c) is the relative horizontal velocity of the flow, 
with reference to the moving frame at speed c. 
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We also recall some basic terminology, as used in [12]: 

• Firstly, the Dubreil-Jacotin transformation maps the unknown domain V to the rect¬ 
angle 

R = ■ -vr < Q < TT, Po < P < 0} , (2.8) 

The independent variables as they appear in the Dubreil-Jacotin transformation m 



—n y=-d X 71 —n P=Pq ^ ^ 


Figure 1: Dubreil-Jacotin transformation 


are 


q = p = —'ip. 


• Secondly, the height function of the wave above B is defined as 

Kq^p) ^y + d. 


(2.9) 


This function satisfies the following nonlinear BVP on the rectangle R. We omit here 
the details of this construction, and refer to m- 

Definition 2.1. The constitutive equations for the height function h{q^p), which is even and 
27v-periodic in q, are the following: 

{ TL[h] := (1 -h hq)hpp - 2hphqhpq + hphqq - = 0 on R, 

Bo[h] := 1 + /i^(g,0) + {2gh - e)/ip(g,0) = 0, (2.10) 

Bi\h] h{q,po) = 0. 

Moreover, the free boundary S is given by the expression 


h{q, 0 ) = p{x) + d. 


In this equation 7 denotes the vorticity as defined in (2.3), which is assumed constant, 
being indicative of underlying currents. Vanishing vorticity is the property typical of uniform 
currents and a constant vorticity characterizes the linearly sheared tidal currents. We also 
recall that g is the gravitational constant of acceleration, Q is the hydraulic head and po is 
the relative mass flux, introduced above. Moreover, we define the height of the wave as the 
maximal variation of the oscillations of the free surface, given by 


a := max h{q,0) — min /i(g, 0 ) =/i( 0 , 0 ) —/i( 7 r, 0 ). ( 2 . 11 ) 

gG[—7r,7r] 9e[—7r,7r] 
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In fact, the BVP (2.10) is a reformulation of the Euler equations restricted to D, under 


the assumptions made previously in this section. 

In the following subsections we discuss some results from m, and provide a new point 
of view, which allows a generalization of these results. 

2.1 Laminar solutions 

The laminar flows are readily obtained as the ^-independent solutions of the boundary-value 


problem (2.10), given by the following formula 

2(p - Po) 




VA - 2'yp + -v/A^-^o^ ’ 
provided that the parameter A > 0 satisfies the equation 

45^01 


Po < P < 0, 


Q — A + 


( 2 . 12 ) 


(2.13) 


\/A + \/X — ‘2j'ypQ 

2.2 Solutions of the linearised problem 

Apart from the laminar flows one can obtain the solution of the linearised problem. For this, 


the BVP (2.10) is linearised around the laminar flow H{p; A). For the speciflc value A* that 
satisfies the dispersion equation 

A , , / 2po 


+ tanh 


= 0, 


(2.14) 


g-'fVX V\/A + yA^^2^, 

an existence result for the solution of this BVP is presented in m- Moreover, the linearised 
solution is given by h{q,p] h) — H{p] A*) -h b m{q^p)^ where 

2(p-po) 


VA* - 2po7 . 1 
m[q^p) = — , , _ - smh 


cosq. 


(2.15) 


x/A* - 2p7 V \/A* - 2'yp + y/X^ - 2'ypo ^ 

In our investigation we start from this result and we obtain an extension of it by making a 
two-fold interpretation. 

Firstly, we interpret the function 

hiq,p-,b) = H{p-X^)+ b m{q,p), (2.16) 

as a perturbation of a laminar solution of the system ( |2.10[ ), in the following sense: the 
system (2.10) is satisfied up to order 6^, i.e., 

'H[h]{q,p) = 0{b^), Bo[h]{q) = 0{b‘^) and Bi[h]{q) = 0 (2.17) 

and the height of the water wave is of order 6, i.e, 

/i(0,0; b) — /i(7r, 0; 6) = 6 [m(0, 0) — m(7r, 0)] = 0{b). 

Having in mind that the wave height vanishes for laminar flows, we can view the expression 


(2.16) as an approximation of small amplitude water waves. 


Secondly, we can regard this expression as the first order asymptotic expansion of the 


exact solution to the BVP (2.10), which leads to the natural question of finding higher-order 
terms. This point of view is motivated by a bifurcation argument, provided in |12] and 
discussed below. 
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2.3 Bifurcation 

We define the curve that represents the laminar flows 


r={(Q(A),if(p;A)):A>0}, 

with Q and H given by (2.13) and ( |2.12[ ), respectively. 

Restricting our attention to rotational flows with constant vorticity 7 , in m it was proven 
that near this curve, as the parameter A varies, there are generally no genuine waves, except 
at critical values A = A* determined by the dispersion relation (2.14). Near this bifurcating 
laminar flow i7*, we have two solution curves: one laminar solution curve A 1 -^ i7(p;A), 
where A and Q are related by ( 2.13[ ), and one non-laminar solution curve Q h{q^p]Q) 
such that hq ^ 0 unless h = i7*, see Figure]^ 



Figure 2: Bifurcation on the curve of laminar flows. 


2.4 Higher order approximation 

Using the above argument we observe that the variation of the parameter Q (which is the 
hydraulic head of the flow) about the uniquely determined value Q* implies an approximation 
of non-laminar flows, in the following sense. Consider the approximation for the hydraulic 
head of the flow, 

m 

n=l 

and the approximation for the height function h{q^p] Q), 

m 

h{q,p;Q) « h^^\q,p\h) = '^hn{q,p;Q*)b'^, (2.18) 

n=0 

with hoiq,p] Q*) = H{p, A*). 

Our goal is to determine {hn{q,P', such that the system 

order i.e., 

n[h^^^]{q,p) = 0(6"*+^), Bo[h^^^]{q) = 0{b^+^) and Bi[h^^^]{q) = 0 . (2.19) 

Moreover, the dominant term of the wave height of this approximation of the water wave 
is of order 6 , i.e, 

aim) ^ Q. Q; 6 ) = 6 [hi(0,0; Q*) - hiiir, 0; Q*)] + 0{b‘^). (2.20) 


( 2 . 10 ) is satisfied up to 
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This condition indicates the accuracy of approximation of genuine waves in the following 
sense: The relation ( 2.19[) sho ws that is an approximation of a solution up to order 
6 ^, whereas the relation ( 2 . 20 ) shows that differs from the laminar solution at order b. 


Thus is ‘closer’ to a non-laminar solution. 

In Section 3, we present the results for m = 2. In Section 4, where we work out the case 
m = 3, the situation is different. There we observe that there is no h^^\q^p]h) such that 
both conditions = (9(6^) and Bo[h^^"^]{q) = (9(6^) are satisfied. This means that 

the approximation made on ( 2.18| ) is accurate up to m = 2 , and there appears an additional 
error for m > 3. Thus, we introduce an auxiliary parameter which allows the control of this 
additional error. Indeed, we construct the following estimates 




< eib^ and 




<€2b^ 


for some small ei > 0 and 62 > 0 , which depend on this auxiliary parameter. 

Moreover, for particular values of the auxiliary parameter we can make each of 61 and 62 
arbitrarily small (but not simultaneously), i.e., 

• either, 'H[h^^^]{q,p) = (9(6^), 

. or, Bo[h^^^]{q) ^ 0{b^). 

The selection of the proper value of the auxiliary parameter is discussed in the last parts of 
Section 4.1 and Section 4.2, as well as in Section 5.2. 


3 Second Order Asymptotics 

The problem of the extension of the asymptotic expansion of the solution of ( 2.10| ) takes 
the following form: with defined by (2.14), i7(p; given by (2.12) and m{q^p) given by 
(2.15), determine u{q^p) such that the function 

h*iq,p-,b) ^ H{p-, X,,) + b m{q,p) + b‘^u{q,p), (3.1) 

satisfies the system ( 2 . 10 | ) up to the second order, i.e. 

T^[h*]iq,p) = 0{b^), BQ[h*]{q) = O(b^) and Bi[h*]{q) = 0. (3.2) 


Note that, in the absence of stagnation points, the real-analyticity of the wave profile for flows 
of constant vorticity, established in uni, ensures the existence of a power series expansion. 
Since the first-order (linear) approximation of the wave profile fails to capture the effects of 
the underlying current, the issue of obtaining high-order approximations should be explored 
in this regard. 


3.1 Irrotational flow 

For the purpose of comparison we analyse first the irrotational case 7 = 0 . This simpler case 
provides us with the opportunity to explain in more detail the strategy we use throughout 
this paper to obtain higher order terms of the expansion of the solution. 




























In the irrotational case, the laminar solution H{p] simplifies to 


H{p-,K) = 


P-Po 
y/K ’ 


the linearised solution m{q,p) becomes 


m{q^p) = sinh ( ^ j cosg, 
where A* > 0 satisfies, in this setting, the dispersion relation 

A + ^ tanh ~ ^ * 


Moreover the corresponding value of Q is 


Q* = A* - 




(3.3) 


(3.4) 


(3.5) 


(3.6) 


We first apply the operator defined in ( |3.2[ ), to the expression (3.1). The outcome 
will be a polynomial in the b variable, of order up to b^. The specific above choices yield 
the vanishing of the zero and first order of the polynomial; for a detailed discussion of this 
aspect, we refer to the HU. Our present goal is the vanishing of the b^ coefficient of the 
polynomial. 

In particular, using the fact that Hq = 0, we get 

T-L[h ] — Hpp + b(jTlpp + HpTTlqq^ 

T b ^Upp + HpUqq Hpp ^HpipfJlqTflqp T 0(5 ). 


Observing that (3.3) yields 




and Hpp = 0 for all p G [po, 0] 


and snbstituting(|3.4[) in the previous expression for Hlh*], we get 


H[h*] = 


Upp + -^Ugg - ^ sinh ( 2 ^ 


+ 0(b^). 


(3.7) 


(3.8) 


Proceeding similarly for the boundary condition, i.e. for the operator in (3.2), yields 

Bom = 1 + (2p77 - Q)Hl + b {2Hp[{2gH - Q)mp + gHpm] } 

+ b^ [rUg + AgmHpmp + {2gH - Q){m^ + 2HpUp) + 2gHpu\ 

+om , 


expression and using the fact that (3.5) is equivalent to the identities 


evaluated at p = 0. Applying the equalities (|3.3|), (|3.4|), (|3.6|) and (|3.7|) in the previous 

and sinh ( ) =- , (3.9) 


cosh 
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we get that 


Bo[h*] = -b‘^2^/X~A {up + 


3/2 


U 


X 


p=0 


- A2 


35 ^ + + ( 35 ^ - A^) cos(2g) 


(3.10) 


+ Oib^). 


Consequently, we obtain the constitutive equations for the function u{q,p)^ even and 27t 
periodic in q variable, 


Upp + -^Uqq = ^ sinh ( 2^—^ 

Ajf; y Y 
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1 


1 


. 3/2 


u = 


4 \/Q 

= 0 , 


2 _ \2 
A>(; 


35 ^ + A^ + ( 35 ^ - A^) cos(2g) 


{q,p) e i?, 


p = 0, 


P = Po- 


(3.11) 


In order to solve the above BVP we use separation of variables and we take in consideration 
the following facts: 

1 f P ~ Po 

• A special solution of the PDE is given by - sinh ( 2 —-=^ 

4 V V Ax. 

• The function u{q^p) is even and 27r periodic. 

• We must have ?i(g,po) = 0. 


Therefore, we get the general form 


u{q,p) = ^sinh + Ao{p - po)+ 


cos{nq) sinh 


n{p-po) 
'S/ Ax; 


(3.12) 


and we have to determine the constants Aq and such that the Robin Boundary 

Condition (BC) at p = 0 is satisfied. Indeed, the separation of variables says that the 
general solution of the Partial Differential Equation 

^pp ~ ^ 

Ax; 

on the rectangle is a bilinear combination of trigonometric functions in q and exponential 
(hyperbolic) functions in p i.e. {cos(ng), sin(ng)}, n G N and exp(±np/vX^), n G N. This 
is due to the periodicity in q^ The evenness in q excludes the dependence on sin(ng), n G N 

and the condition = 0, on p = po gives exactly the sinh dependence. 

Thus, applying the form (3.12) to the second of the equations (3.11), we get 

Ax. 3p^ - A^ 1 


An = 


4 (^2 _ \2 . 3 / 2 ’ 


A 2 = = 0, n > 3, 


with Ax. given as the solution of (3.5). 


^Once we have obtained {cos(ng), sin(nq)}, separation of variables gives the exponentials exp(±np/vV). 
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3.2 Constant Vorticity 

The strategy is conceptually the same for the rotational case when the vorticity 7 is constant. 
The first goal is to derive the analogue of the BVP (3.11). In this case, the laminar solution 
H{p) is given by 

p-po 


where 


= 2 - , 
r{p) + r{po) 

r{p) = a/A* - 2jp, 


(3.13) 


(3.14) 


with A* > 0 satisfying the dispersion relation ( |2.14[ ). Moreover, the linearised solution m{q^p) 

(3.15) 


is given by (2.15) and rewritten as follows 

ripo) 


Tn(q^p) = —sinh (H(p)) cosq. 
r{p) 

The corresponding value of Q, using ( |2.13[ ), is 


Q* = A* - 


y/X^ + \/A>t, — 27P0 


Then we evaluate the system (2.10) for /i*, defined by 

h*{q,p;b) = H{p) + b m{q,p) + h^u{q, p ), 


(3.16) 


(3.17) 


with H{p) and m{q^p) given by ( |3.13[ ) and (3.15), respectively. Demanding the validation of 
this system up to the second order, i.e. 


n[h*]{q,p) = 0{b ^), Bo[h*]{q) = 0{h^) and Bi[h*]{q) = 0, 


(3.18) 


yields the analogue of the BVP (3.11). After tedious but straightforward calculations we 
get the following form for the BVP, for the even and 27r-periodic function u{q^p) on the 
rectangle R = [— 7 r, 7 r] x [po^O]: 

1 37 

Upp{q,p) + -^^Uqq{q,p) - -j^Up{q,p) = fo{p) + f 2 {p) cos( 2 g), 


r^{p) 


r^{p) 


9 


Upiq, 0 )- ^u{q, 0) = go + 52 cos( 2 g), 

A* 

u{q,Po) = 0, 


(3.19) 


where 


f 1 A r'^iPo) 

/«W= ,7(p) 

37 = 37 ( 7 ^+ 2 r 2 (p)) cosh 


+ 2r(„) ( 37 ^ + 2rHp)) sinh ( 2 "-’'<»>) 1 

•yr^ipo) 

[ 37 ^ 2r\p) (37^ + 4r2(p)) cosh ~ 


+ 67 r(p)sinh( 2 ’'<*’>-’'<«'>)], 
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r^jpo) 3g^ + Xl 

2A* -\ l - ' ys/K (25 - 7\/^) ’ 

r^(po) 3g^-Ag 

2A* g"^ -\l- 7v^ (25 - 7v^)' 


Remark 1 . The above BVP (3.19) reduces to ( |3.11[ ) for 7 = 0 , 
vanishes when 7 vanishes. 


i.e. 


the funetion f 2 {p) 


Applying separation of variables to (3.19), we obtain the following formulation of the 
above BVP: Determine u{q^p) given by 


u{q,p) = Uo{p) + E- n{p) cos{nq), 


(3.20) 


n=2 


where the functions Un{p) satisfy the following BVPs: 


n- , . 37 , 


^nip) - fT^^nip) - ^ fnip), P^ (PO, 0) 


<(0) 


3/2 


a: 


^n(O) — Pn, 


(3.21) 


UniPo) = 0 , 


where fo{p), f2ip), go and g 2 are defined above. Moreover, fn = 0 and gn — 0, for n > 3. 
Now, we make the following change of variables 


Unip) = 


Unirip)) 


with rip) — a/A* + 2'yp 


rip) 

which transforms the above BVPs into the following: 

72 ^2 

r e (ro, VV 
r = a/V 

r = ro := r(po) = ^A* + 27^0- 


(3.22) 


(f 71“ 

-^Unir) - ^Unir) = fnir), 

d , 5 + 7^V, 

-—1X72 - Upl — ppij 

ClT 7 


(3.23) 


'^n — O 5 

Here fn and pn are defined by 

fnir) = fnir ip)) = ^^^fnip) 
T 


and 


/V 

9n — Qn-) 

7 


n G N. 


The solution basis of the homogeneous differential equations appearing in BVP (3.23) is 
given by 


{ 1 , r} for n = 0 and 


. <!sinh(n^^^bcoshfn'’ 


7 


7 


for n > 1 . 
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Using the method of variation of parameters, we obtain the special solutions of the BVP 


(3.23). In particular, we seek for special solutions of the inhomogeneous problems of the 
form 

• For n = 0: 


under the constraint 


«oW + I^Q{r)r = 0. 


• For n > 1: 


u^\r) — an{r) sinh [ n-—— ) + /3n(^) cosh I 


7 


7 


under the constraint 


7 


sinh ( n-—— ] + /3^(r) cosh [ n-—— ) = 0. 


7 


The insertion of these forms in the inhomogeneous differential equations appearing in BVP 
(3.23) yields {an^Pn}^ n G N. Consequently, the special solutions read 


= -7?’^ 


1 — cosh 2 


J - rp 
7 


sinh I 2 


8r^ 


- rn 


r - rp 
7 


and 


1 — cosh 2 

_ .,^2 V 7 




r — tq 


sinh 2 


8r2 ^0" 

The next step is to apply the linear combination 


4r 

r - ro 


7 


4r 


Un = 


to the Robin condition of (3.23) and obtain the specific value for E^. One can readily see 
that = 0 for n > 3. 

Finally, after determining Eq and £^2, we use (3.22) and, by means of (3.20), we obtain 
the following formula for the second order term of the expansion (3.17): 

u{q,p) = r‘^{po) 


H{p) 1 — cosh {2H(p)) sinh(2i?(p)) 

^O-rTTT “ 7- ^ Q/ X -1- 


r{p) 


8r^{p) 


4^2 


+ r^(po)cos(2g) 


sinh {2H{p)) 1 — cosh {2H{p)) sinh {2H{p)) 

^2-TX-7- ^ Q/ N -K 


where 


Co = 


r{p) 8r^{p) 4r^(_p) 

[^9 {9 ~ 7V^) + (7^ “ ^*)] (^(Po) + V^) 


(3.24) 


{g-jVK) -A2 [2gpo +X^r{po) + ^/Kr‘^{po)] 
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and 


C2 = 


‘^9 ~ 7 a/^) + ( 7^ ~ 3Ax,) 


8A 


5/2 


Therefore, the function /i* given by (3.17) satisfies the system (3.18). We note that for 


the value 7 = 0 , the expression (3.24) takes the form (3.12), meaning that we obtain the 


irrotational case as a special case of the constant vorticity formulation. 

4 Third Order 

With a similar methodology as that used to derive approximations of second order, one can 
derive approximation of third order. The results are provided below; the reader can easily 
verify them with some Computer Algebra Software. 

In particular, we use the formula 


b) = H{p) + b m{q,p) + b^u{q,p) + b^w{q,p), 


(4.1) 


with i7, m and u given in the previous section and A* given by (2.14). Then our aim is to 


determine w so that h satisfies the system (2.10) up to the third order, i.e. 

'H^]{q,p) = 0{b ^), BQ[h]{q) = 0{b‘^) and Bi^]{q) = 0. 


(4.2) 


However, there is no such h{q^p]b). Thus, instead, we introduce an auxiliary parameter B 
(for the irrotational case we denote it by B) which allows the control of this additional error 
estimate. Indeed, we construct the following estimates 




< eib^ and 


Bo[h^^^]iq) 


< e2b\ 


for some small 6 i > 0 and 62 > 0 , which depend on this auxiliary parameter. 


4.1 Irrotational case 

Following the organisation of the previous section we first discuss the case that 7 = 0 . 
Therefore, if we apply ( |4.1[ ) to the system (2.10), we get a polynomial on b of finite order (in 
particular of order b^). The requirement that the coefficient of the b^ term vanishes yields 
the following BVP for the even and 27r-periodic (in the q variable) function w{q^p): 


Wpp + ^Wqq = fi{p) cos q + /a(p) cos(3g), 

A>k 


W, 


'p > 3/2 

A>k 


w = gicosq + gz cos(3g), 
re = 0, 


{q,p) e R, 

p = 0, 
P = Po, 


(4.3) 


where fi{p) and / 3 (p) are known functions (bilinear combinations of hyperbolic functions and 
polynomials) and gi and are known constants, dependent on the parameters p, po and A*. 
We use the fact that this problem has the same general solution (for the homogeneous case) 
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with the problem (3.11) and we determine its special solutions by the method of variation 
of parameters. Then the expansion 

Hq,P; b) = H{p) + b m{q,p) + h^u{q,p) + b^w{q,p), 

with w{q^p) given the following formula 


w{q,p) = 

+ 


Bi {p — po) cosh ( ^ ) + -^2 sinh { 3^ 


■\/V 


Di sinh ( ] + D 2 sinh { 3^ 


with 




A* 


4 o2 _ \2 ,3/2’ 

^9 qpo + XJ 




B2^ 


vV 

cos(3g), 


cos q 


(4.4) 


32Ar 


D,=- 


3g^ - A^ 
32A2 


and = + 


32A2 


32A2 


satisfies the conditions 


'bi[^{9,p) — 0{b‘^) and Bi[h]{q) — 0. 
For the Robin boundary condition we take 

= -Bocosq b^ + 0{b'^), 

with the constant Bq given explicitly by the expression 


(4.5) 

(4.6) 


Ba = 


1 


1 


Vpo - A^)^ g{3g‘^ + xl) 


2^2' 


I \3/2 

gpo + A*' 


2A* (^g2 _ a 2)3/2 
Thus, we propose the following formulation 

w{q,p) = w{q,p) - B{p-po)cosq, 


2A* 


(4.7) 


for some constant B. Then the condition ( |4.2| ) changes accordingly. In particular, ( |4.5| ) 
becomes 

'H[h]{q,p) — B ^ cosq b^ + 0{b‘^) (4.8) 


A,); 

and the Robin boundary condition takes the following form: 

2 


Bo[h]iq) = 


^ [gPo + A*/") B-Bo 


cosq b-^ + 0{b^ 


(4.9) 


Therefore we can choose some appropriate value of B in order to minimize the errors observed 
in the above conditions. For example, 

• S = 0 ensures T-L^] = (9(6^). 

A* 1 


• B — Bci 


gpo + A: 


3/2 


ensures B{)[h] — 0{h^ 
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We observe the magnitude of the error of the term, for the operators H and Bq, to be 

(4.10) 


Sd := > 0 


and 


respectively, with 


£c So ~ ^ {qP^ + ^ S > 0, 

A* 1 


(4.11) 


0 < S < So- 


I \3/2‘ 
gpo + Axe 


4.2 Constant vorticity 

While the situation is similar for this case, it is, however, much more complicated. Indeed, 
we apply the expansion of the form (4.1), with i7, m and u given by (3.13), (3.15) and (3.24), 


respectively, to the system ( |4.2| ) and we get the following BVP on the rectangular i?, for the 
even and periodic function w{q^p): 


1 37 

Wppiq,p) + :^j^Wqq{q,p) - ^^j^Wp{q,p) = fi{p) cos q + fsip) cos(3g), 

Wpiq, 0)-0) = 9i cos q + gs cos(3g), 

Ax; 


(4.12) 


U!{q,Po) = 0 , 

where 

• flip) 2 tnd fsip) are known functions explicitly dependent on r(p), 

• gi and are known constants explicitly dependent on the parameters 7 , Ax=, po and 

for which we avoid to write the complete expressions for matters of brevity. Proceeding 
as in Section 3.2, we reduce the solution of the above BVP to the following formulation: 
Determine re(g,p), of the form 


n=0 


with Wnip) given by 


Wnip) = 


w{q,p) = ^Wn{p) cos(ng), 

Wn{r{p)) 
r(p) 


with r(p) = i/A* — 2jp, 


(4.13) 


(4.14) 


and where Wn(r) are the solutions of the BVPs 

7i?‘ 

- ::^Wn{r) = fnir), 
d ^ Q- 7 W , _ . 

Wfi + — g^^ 

dv 7 Ax; 

Wn = 0, 

for some known functions fnif) and constants Qn- 


r G 


(rQ,V^* 

r = a/V 
= ro = y/A* - 27P0, 


(4.15) 
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Remark 2. The functions /n(^) are bilinear combinations of hyperbolic functions and poly¬ 
nomials in 1/r. The constants arc dependent on 7 ^ g, and tq. 


Eventually, using the arguments of the previous subsection we derive the following for¬ 
mula for w{q^p): 


/ \ \ r, 2 ^ \H(p) 

w{q,p) =w{q,p) - Br {po)—^cosq, 

r{p) 

for some constant B and w{q^p) given by 

sinh (i?(p)) + A 2 {p) sinh (3i7(p)) 


-\-Bi{p) cosh {H{p)) -h B 2 {p) cosh {3H{p)) 


cosq 


+ 


A^{p) sinh (^H{p)) -\- A 4 {p) sinh (3H{p)) 

+Bs{p) cosh {H{p)) -h B 4 {p) cosh (3i7(p))j cos(3g), 


(4.16) 


(4.17) 


where {Aj{p)^ Bj{p)}^^-^ are functions of the form 


E Unr" ip) and ^ ^ {p ), 


r^{p) 


n=0 


n=0 


respectively. Here 


^ 2 s 

«i,o = -^7ro 


on,i = 0, 

01,2 = - j^o (35'(5f - 7VV - A*) VV {g - 7W + A*) , 

+?’o ( - - 155(^7 W - 57 (7^ - A*) A*^^ + 3 gfA* (57^ ~^ '^*) )) }/ 

|32 (—ro [g — 7 \/A 7 ) + gV^) {g^ ~ ‘^giV^ + (7^ “ A*) A*) |, 

^0 (35^ - ^giVK + (7^ - 3 A=^) A*) 

, ^1,4 6 , 

7^0 {^ 9 ^ ~ 3 ^ 77 ^ + (7^ - 3 A*) A*) 


Ol,3 = — 
^1,0 = 0, 
^>1,3 = 
biA — 


biA = - 


32A 

97^0 

32 ’ 


5/2 


&1,2 = — 


32A^/^ 


rQ\/V {—‘ig^ + 35'7\/V + A* (— 7 ^ + A=^)) 

4 {gro -{g + iro) v^) (5^ - “^giVK + (7^ - A*) A*) ’ 
(35f^ - Sg'jVK + ( 7 ^ - A*) A*) 


4 {gro - {g + 7 ^ 0 ) \/V) {g - iVK - A*) [g - j^/K + A*) ’ 
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«2,3 = 


t>2,2 — 


03,3 = 


04,3 = 


^4,2 = 


3 2 3 

9/^0 

tt 2 ,l = 0 , a 2,2 = , 

CO 

000 

to 

1 

35 f 7 V% + (7^ - 



32 a!''= 

5 

0 , ^> 2,1 

9 3 

= ’ 


7 ?’o (35^ - 

Sg'jy/K + (7^ - 

3 A^) 


32 Ay^ 

5 

3 2 3 

^ 3,1 = 0 , as ^2 

_ ^0 

32 ’ 

1 

(M 

CO 

coo 

1 

- 3 ^ 7 ^^ + (7^ - 

- 3 A>t;) 


32 A*/^ 

5 

0 , (> 3,1 

37 r^ 

32 ’ 


1 

ooo 

^co" 

to 

- ^gi'/K + (7^ 

— 3 A^) A^) 


32 A=>" 

? 

1 2 3 

3 rQ 

04,1 = 0 , a 4^2 = , 

CO 

ooo 

to 

1 

3 g^^/K + (7^ - 

3 A>(;) A>t;) 

1 

coo 

Z 2 xf- 

5 

279 * 7 ^ + II9 

^ (37^ — 2 A* 

+5A* (77^ - 357 ^A* + 13 A^) - 7 A*^^ 

{64 (g - 

-I'/K) A^l, 


0 , (>4,1 

3 3 

= ’ 


7 >’o (35^ - 

“^giVK + (7^ - 

3 Ax;) X^^ 


02,4 = 0 , 


t>2,3 — ^ 2,4 = 0 , 


03,4 = 0 , 


^ 3,3 = ^ 3,4 = 0 , 


,3/2 


32A 


5/2 


^ 4,3 = ^ 4,4 = 0 , 


where we used the abbreviation r{po) — rg. Consequently, the conditions (4.2) take on the 
following form: 

(4.18) 


'H[h]{q,p) = Br^^ipo)^^ cosq + 0{b‘^) 


and 




2 / N “^gpo + r{po){r{po) + v/A)() 


2r {po)- 


with 


(r(po) + V^) 

X)n=o “Ax 7" 


B-Bc 


cos q b^ + 0(b^), (4.19) 


Bo^n 


- A2 - 7V^ (2g - 7\/X) ’ 
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where 


1 


and 


D = - 


“^gPo + r 


rjpo) + \/A^ _ 

{po)VK {r{po) + ^/K) 


Co = 3gy/K {9g‘^ - Mg'^Xl + 5A^) - 4g^/X^ro {9g'^ - 9g‘^\l + 4A^) 

- ro ~ ~ ’ 

Cl = -A*^^ {2>‘4g'^ - S9g^\l + 4At) - ro (-GO^^A.^ + dSg'^A^ - 7A®) , 
C2 = -3 (^-7g^\‘l + 55rA^ + 2gro\l^^ {9g'^ - 4A^)^ , 

C3 = -A*^^ {Tg"^ - 2A2) - ro {-28g'^Xl + 5A^) , 


C 4 = -SgroXl^'^ + gX^, cs = roA^ . 

Similarly to the previous section we observe that the magnitude of the error of the term, 
for the operators H and Bq^ to be 


£d := > 0 

r{po) 


and 


Sc := Bo - 2r (po) 


2 / 1 ‘^gpo + r(po) W (r(po) + W) 


xl^^ (r(po) + W) 


B>0, 


respectively, with 


(4.20) 


(4.21) 


0 <B <Bq 


)?J‘^ _ r{po) + \/V _ 

2 r 2 (po) 2gpo + r(po)^V (r(po) + W) 


Remark 3. One can see that the case 7 = 0, equivalent to r{p) — x/X^ for all p G [po? 0], in 
the above expressions yields the formulae of the previous subsection. 


5 Illustration of the solutions 

Considering the figures depicted in this section, we recall that the value of the gravitational 
constant of acceleration is p = 9.8 and we fix the relative mass flux po = —2. This section 
is organised as follows. In the first subsection we illustrate the results for the second order 
asymptotic expansion. In the second subsection we present the analogue results for the third 
order expansions. In the third subsection we present some figures that compare the obtained 
second order approximation results against existing results i.e. first order asymptotic expan¬ 
sions. Finally, in the fourth subsection we compare the results between the second and the 
third order approximations. 

The illustrations below are given for several values of the vorticity 7 , as indicated in each 
figure. Moreover, each streamline is defined by [p = constant], for constants ranging from po 
to 0. We note that for the wave height h and the velocity field v^ we have that: 
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^ P — Po gives /i = 0 and 'u = 0. 

• /i is an increasing function with respect to p. 

• 'L’ is an increasing function with respect to p, for g > 0, and a decreasing function with 
respect to p, for q < 0. 


whereas, for the pressure P: 


• P{q,0) = Patm{q) = Patm- 


• P is a decreasing function with respect to p. 


Let us recall the definition of the wave height as the maximal oscillation of the free 
boundary, /i(0,0) — /i(7r,0), with the understanding the the wave crest is located at g = 0. 
Moreover, we recall from m that: 


• The velocity field is given by (c — v) 



• The water pressure is given by 


1 + Q 

P = Patm - - gh--fP+^, 


for some constant Q. 


5.1 Second order 


In this subsection we depict the results that we obtain from the second order asymptotic 
expansion of the solution /i(^,p), given in (3.17). Concerning the other parameters, we 
remind that A* varies with the vorticity according to the dispersion relation (2.14). 

Firstly, for each value of the constant vorticity 7, we have to fix the value of b by allowing 
a small error on condition ( |3.18| ), i.e. 


max{||P||2, ||>So|l2} — ^ for some small e > 0. 
b 

In Figure we show how b varies for different values of vorticity. The value b is increasing 
when 7 G [—4.5,—1.5], with ||P||2 dominant in the above definition. On the contrary, the 
value of b is decreasing when 7 G [—1.5,3], with IIP0II2 dominant in the above definition. 
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Figure 3: The value of b for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 


5.1.1 Wave profiles 




Figure 4: The free surface for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 
Here the wave profile r]{x) is illustrated, which is oscillating around y = 0. 
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Figure 5: The mean depth for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 



Figure 6 : The height of the water along streamlines for 7 = —1.5, over two wavelengths. 
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Figure 7: The wave height for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 
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5.1.2 Velocities 




Figure 8 : The vertical velocity v on free surface for different values of vorticity, varying from 
7 = —4.5 to 7 = 3. 



Figure 9: The vertical velocity v of the flow along streamlines for 7 = —1.5. 
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Figure 10: The horizontal velocity c — u on the straight line below the crest, i.e. q — {) and 
P ^ [:P 05 0], for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 


5.1.3 Pressure 




Figure 11: The water pressure on the bottom, p — p^^ for different values of vorticity, varying 
from 7 = —4.5 to 7 = 3. 
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0.0 

Figure 12 : The water pressure of the flow for 7 = —1.5: {P — Patm) vanishes on the free 
boundary p = 0 and takes the maximal value on the bottom p = po- 


5.2 Third order 

Here, we proceed as follows: 117^112 and ||)Bo|l 2 for the third order asymptotic expansion are, 
in general, functions of 7 , b and H, see S and ( |4.17[ ). So, for each value of 7 , we pick the 

f ecific value of b that we have obtained in the previous subsection, as illustrated in Figure 
Then, for each pair ( 7 , 6 ), we determine the value of B that minimizes the quantity 


11^112 +Poll: 



Figure 13: The value of B for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 


In Figure 13 we show how B varies for different values of vorticity. In what follows we 


will display the analogue of the figures of the previous subsection, but now for the third 
order asymptotic expansion 
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5.2.1 Wave profiles 



- 0.10 



Figure 14: The free boundary for different values of vorticity, varying from 7 = —4.5 to 
7 = 3 . 
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Figure 15: The mean depth for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 
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Figure 16: Each of the disjoint curves displays the wave profile r]{x)^ for different values of 
vorticity, varying from 7 = —4.5 to 7 = 3. 



Figure 17: Each of the disjoint curves displays the free boundary h{q^0) = r]{x) + d, for 
different values of vorticity, varying from 7 = —4.5 to 7 = 3. This illustration, compared to 
Fig. (16), takes into account the mean depth d, which varies with the vorticity 7 , see Fig. 
151 
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Figure 18: The height of the water along streamlines for 7 = —1.5, over two wavelengths. 
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Figure 19: The wave height for different values of the vorticity, varying from 7 = —4.5 to 
7 = 3 . 
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5.2.2 Velocities 




Figure 20: The vertical velocity v on free surface for different values of vorticity, varying 
from 7 = —4.5 to 7 = 3. 



Figure 21: The vertical velocity v of the flow along streamlines for 7 = —1.5. 
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Figure 22: The horizontal velocity c — u on the straight line below the crest, i.e. q — {) and 
P ^ [:P 05 0], for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 


5.2.3 Pressure 
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Figure 23: The water pressure on the bottom, p = po? foi* different values of vorticity, varying 
from 7 = —4.5 to 7 = 3. 



Figure 24: The water pressure of the ffow for 7 = —1.5: {P — Patm) vanishes on the free 
boundary p = 0 and takes the maximal value on the bottom p = po- 
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5.3 Comparison with existing results 

Here we compare the results that we obtained from the second order asymptotic expansion 
with the ones already existing in the literature, i.e. the first order asymptotic expansion. In 
order to choose the value of b for the depiction of the first order asymptotic expansion, we 
proceed as in Section 5.1 and we obtain qualitative results similar to Figure but now b is 
in general smaller, as depicted in Figure 



Figure 25: The value of b for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 

This choice of b indicates the comparison on the wave height for different values of 
vorticity, depicted in Figure There we observe that for large values of the absolute value 
of vorticity, jyj, the wave height is relatively small and the first order approximation shows 
similar results with the second order one. However, the shape of the water profile differs 
significantly, see Fig. and Moreover, for the rest values of vorticity, we observe a 
discrepancy in the wave height between the two approximations. 
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Figure 26: The height of wave for different values of vorticity, varying from 7 = —4.5 to 
7 = 3. 


5.3.1 Negative vorticity 7 = —1.5 

In what follows we fix our vorticity to 7 = —1.5 in order to illustrate the improvement 
provided by higher-order asymptotics. 



First order 

Second order 


Figure 27: The free surface for vorticity 7 = —1.5. The mean depth is 0.80072 and 0.79222 
for the first and the second order approximation, respectively. 




(b) Second order asymptotic. 


Figure 28: The height of the water along streamlines for 7 = —1.5, over two wavelengths. 
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First order 

Second order 


Figure 29: The vertical velocity v on the free surface, for vorticity 7 = —1.5. 



First order 

Second order 


Figure 30: The water pressure: P — Patm on fhe bottom, for the choice of vorticity 7 = —1.5. 


5.3.2 Positive vorticity 7 = 1.5 

In what follows we fix our vorticity to 7 = 1.5 in order to illustrate the effect of higher-order 
asymptotics. 



First order 

Second order 


Figure 31: The free surface for vorticity 7 = 1.5. The mean depth is 0.758042 and 0.752228 
for the first and the second order approximation, respectively. 
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(a) First order asymptotic. (b) Second order asymptotic. 

Figure 32: The height of the water along streamlines for 7 = 1.5, over two wavelengths. 



Figure 33: The vertical velocity v on the free surface, for vorticity 7 = 1.5. 



First order 

Second order 


Figure 34: The water pressure: (P — Patm) on fhe bottom, for vorticity 7 = 1.5. 


5.4 Comparison of second and third order asymptotics 

In what follows we fix our vorticity to 7 = 1.5 and we compare the difference between the 
second and third order expansion. 


In Figure 36 we depict, for different values of vorticity, the curves hqq = 0 for p G [po, 0], 
thus indicating where the curvature changes. Indeed, let us fix the value of vorticity 7 . If 
for instance p = p is fixed, then the equation 


d^h{q;p) 

dq^ 


= 0 


has a unique solution q — q, where q G ( 0 , 7 r/ 2 ). Allowing p varying from po to 0, we get the 
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Figure 35: The wave height for different values of vorticity, varying from 7 = —4.5 to 7 = 3. 

curve p q{p). For different values of the vorticity we get the different curves depicted in 

Figure observing where the wave profile changes from convex to concave. 




(a) Second order asymptotic. (b) Third order asymptotic. 

Figure 36: The curve tracking the change of curvature, for values of the vorticity varying 
from 7 = —4.5 to 7 = 3. 


6 Conclusions 

Among the conclusions that we can extract from our asymptotic expansions and the associ¬ 
ated figures, the most important are the following: 

• We observe a considerable improvement from the first to the second asymptotic ex¬ 
pansion, noticeable in the depictions of the wave height in Fig. as well as in the 
depictions of the profile of the wave in Fig. [^and[3^ 
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• We observe a small difference in wave height from the second to the third asymptotic 
expansion, cf. Fig. 1^ However, there is a considerable difference in the change of 


the curvature, see Fig. 36 This is due to the choice of the parameters b and B. In 


the particular examples we did choose to keep the same values for b and to minimize 
the errors 117^112 and ||)Bo||2 5 meaning that we kept approximately the same wave height 
and improved the accuracy of the wave profile. 


• Fig. [T^ shows that an opposing current has a steepening effect on the wave profile. As 
already pointed out in the Introduction, first-order linear theory does not capture this 
important feature of wave-current interactions. 


Asymptotic expansions beyond those obtained here will be presented in upcoming pub¬ 
lications. 
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